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Using the data sample of 10.5fb 1 collected by the Belle detector, we searched for two-body charmless decays involving vector 
mesons. The clear signal signature of B^ — » ifiK^ is seen and its branching fraction is measured to be (1.39^'gg ± 0.25) X 10~ 5 . The 
evidence for — > p°7r , B° — > p^^i an( j fl± tfcK*^ 1 is reported and their branching fractions are determined. No significant 
signals are observed for B° -> -> -> pT^± i fl± _> p°K ± , - 



and — > ojX^ only the 90% C.L. upper limits are 



given. 

1 Introduction 

The study of charmless B decays will play a significant 
role in checking the KM matrix description and in test- 
ing the validity of the Standard Model description of CP 
violation. 

Decays involving b— > sss transitions such as B— > ipK 
cannot occur via a tree process and are expected to be 
dominated by the penguin diagram, thus their measure- 
ment will give information on the strength of the penguin 
transition. They are also sensitive to physics beyond the 
Standard Model, since non-SM particles can contribute 
via additional loop diagrams. 

The decay B° — > p ± 7r T is one of the preferred modes 
for determining the angle (f>2 of the unitarity triangle. 



The decays B d 
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and B d 



cow are interesting as 



they may exhibit large direct CP violation. The B ± ^> 
p°K ± channel may yield information about the angle (j>s . 

The first searches for B — > 4>K decays modes were 
performed by the CLEO collaboration in data samples 
of gradually increasing integrated luminosities of 2.42, 
3.11 arid 5.5 fb _1 but no evidence for these modes was 
found p. The first evidence for the decay J5 ± — » (pK^ 
in a data sample of 5.5 x 10 6 BB pairs was reported by 
the Bella collaboration at the ICHEP2000 Conference 
in Osakaa. The branching fraction was measured to be 
B(B ± ^(f>K ± ) = {VT2toil ± 0.18) x 1(T 5 . At the same 
conference, CLEO0 presented its first hints of the obser- 
vation of this decay mode, quoting a branching fraction 
of (0.64±°;25+°'^) x 1(T 5 . Recently CLEO reported! 

measurements based on the data sample of 9.7 x 1Q 6 BB 
pairs: B(B ± -> c^K^) = (0.55±J^ ± 0.06) x 10~ 5 and 

B(B ^ <f>K* a ) = (1.15lg;g^"J) x 10- 5 and 90% C.L. 
upper limits for the branching fractions of J f?* ± — ^ (jjK*^ 
and B° — > (hK° decays of 2.25 x 10~ 5 and 1.23 x 10" 5 
respectively □. 

Charmless B decays into modes with p(ir/K) and 
uj(tt/K) are hampered by large combinatorial back- 
grounds. These modes have only been seen recently, first 



Table 1: The measured branching fractions by other experiments 
[ID" 5 ]. 



Mode 


BABAR(BF) 


CLEO(BF) 


B±^ p ^ 


2.4 ±0.8 ±0.3 


1.04 ±0.33 ±0.21 


B±^ pOR* 


1.0 ±0.6 ±0.2 


0.84 ±0.38 ±0.18 




4.9 ± l-3±g;f 


2.76±g;|t ± 0.42 


B ^ pTK± 




< 3.3 (90%C.L.) 


B±^ W7T ± 




1.13 ± 1.4 






< 0.79 (90%GL.) 


B±^ uuh ± 


0.89 ±0.54 ±0.22 





by CLEOB and then by BAB ARB. The current branch- 
ing fractions measured by CLEO and BABAR are given 
in Table |. 



2 The data sample and the Belle detector 

The presented results are based on the data set col- 
lected by the Belle detector 13 at KEKBEl, the asymmet- 
ric B-factory at KEK. It consists of 10.5fb _1 taken at 
the T(45) resonance and 0.6fb _1 taken below the BB 
production threshold used for continuum studies. The 
Belle detector is a general purpose magnetic spectrome- 
ter equipped with a 1.5 T superconducting solenoid mag- 
net. Charged tracks are reconstructed in a 50 layer Cen- 
tral Drift Chamber (CDC) and in three concentric layers 
of double sided silicon strip detectors of the Silicon Ver- 
tex Detector (SVD) . The SVD allows precise reconstruc- 
tion of secondary decay vertices. The charged particle 
acceptance of the spectrometer covers the laboratory po- 
lar angles 17° < 8 < 150° which corresponds to ~ 90% 
of the full CMS solid angle. Photons and electrons are 
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identified using the CsI(Tl) Electromagnetic Calorime- 
ter (ECL) located inside the magnet coil. Muons and 
i^£'s are detected using resistive plate chambers embed- 
ded in the iron magnetic flux return (KLM). Charged 
particles are identified using specific ionization losses in 
the CDC and identification information from the Aerogel 
Cherenkov Counters (ACC) and Time of Flight Coun- 
ters (TOF). By these three methods, K/n separation is 
achieved over the momentum range from about 0.2 to 
3.5 GeV. The high momentum kaon identification effi- 
ciency is determined to be ~ 80% (at a purity of greater 
than 85%). 

3 Event selection 

In the search for the above channels we take full advan- 
tage of asymmetric e + e~ collisions resulting in boosted 
BB pairs, good vertexing and particle identification ca- 
pabilities of the Belle detector. 

The details of analysis are optimized for each channel 
separately. For lack of space, only the main steps in the 
analysis can be described. 

We search for <f> — > K + K~ decays by selecting pairs 
of oppositely charged tracks consistent with kaon hy- 
pothesis (which accepts kaons with > 90% efficiency) 
and \Mkk — < lOMeV. Events with a candidate 
4> are accepted if the 4> momentum in the CMS exceeds 
2.0 GeV. 

Candidate p° — > 7r + 7r~ decays are required to have 
dipion invariant masses which satisfy \M nn — M p a\ < 
220MeV, while the p^ candidates are required to satisfy 
\M n ± n o -M p ± \ < 200MeV. 

The uj is reconstructed in decay mode uj —* 7r + 7r~7r°. 
We require both 7r°'s daughter photons carrying energy 
greater than 50McV. A variable, defined as the cross 
product of the 7r + 7r~ momentum vs 7r° in the uj rest frame 
(called uj -amplitude ) is used to suppress the combinato- 
rial background. The 7r + 7r~7r° mass is required to be 
within a ±30MeV window from the nominal uj mass. 

The charged tracks used in the oj/p reconstruction 
are required to be inconsistent with lepton or kaon. 
The required level of particle identification (PID) for ac- 
companying kaons K^, K*^ , K*°, is optimized using 
the measured kaon efficiency and pion misidentification 
curves. At the chosen cut value the kaon identification 
efficiency exceeded 80%. For Kg — > tt + tt~ decays we use 
oppositely charged track pairs where the displacement 
of the 7r + 7r~ vertex from the interaction region in the 
transverse (r — <j>) plane is more than 1 mm. The <f> co- 
ordinate of the 7r + 7r~ vertex point and the <f> direction 
of the 7r + 7r~ vertex point agree within 0.2 radians and 
— m K a\ < 15MeV. The K* candidates are recon- 
structed in four modes: K*° -> , K*° -> K s ir°, 
K *± _^ K ± n o and K *± _^ K^ir*. The invariant mass 
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Figure 1: The M(, (top) and AE (bottom) distribution for selected 
— > 4>K^ candidates. The curves are projections of the un- 
binned ML fit. 

is fitted with the beam spot constraint (described below) 
and the signal window is defined within ±50MeV of the 
nominal resonance mass. 



4 B meson reconstruction 

The selected vector meson and the pseudoscalar (or vec- 
tor) are combined to form the B-meson decay candidate. 
The vertex fit of the B candidate is performed with a 
beam spot constraint (a x w 100 ^m , <r y w 5 pm and 
c z « 3 mm) enlarged by a halo corresponding to the 
width expected for the B meson lifetime (rs 20 pm in 
the transverse plane) . A run dependent interaction point 
position and beam spot size are used. 

The candidate B decay is identified by a beam con- 
strained mass Mf, = \J (y / s/2) 2 — P^ 2 and the calculated 
energy difference AE = E* B - y/s/2, where E* B and P B 2 
are the energy and 3-momentum of the B candidate in 
the T(4S) rest frame. The width of the B signal region 
in Mb and AE depends on the expected resolutions for 
each decay mode. 

The continuum background has been suppressed us- 
ing observables: the angle between the B candidate mo- 
mentum vector and the thrust vector of the remaining 
tracks of the event (| cos(8thr-B)\), and the B candidate 
production angle in the CMS (| cos(^)|). For B -> PV 
decay channels an additional variable is used: | cos(##)|, 
where the helicity angle 9h is the angle between the di- 
rection of the K/tt and the momentum vector of </>/ ' p/ui, 
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Table 3: The fit results for _B° -> o+ttT modes. 
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Figure 2: The M b (top) and AE (bottom) distributions for _B ± — > 
cpK^* candidates. The curves are projections of the unbinned ML 
fit. 



data sample 


fitted signal yield 


significance 


cos(Q H ) > 


10.9±|;1 


2.6a 


cos(Oh) < 


10.1±|;| 


2.4a 



linear function. 

The accepted candidates for B^ — > (pK^ are shown 
in Fig. [jj The fit gives 17.8±l* si S nal events. The sta- 
tistical significance of the signal yield is 7.5o[]. The effi- 
ciency for this decay mode is 11.8%, which includes de- 
tector acceptance and the intermediate BF. 

The results of the fit for B^ — > (pK*^ decays can- 
didates is shown in Fig. || . The fitted number of signal 
events is 6.5j^'7 with significance of 3.6a . The efficiency 
for this decay mode is 4.0%. 

The results of the fits for other B — > cf>K modes are 
summarized in Table |^. 



Table 2: The fit results for B —> 4>K modes. 



Mode 


fitted signal yield 


significance 


B ± -> (j)^ 


17-8±t5 


7.5a 


B° -> cj)K* 


6.5i| : 5 7 


3.6a 


B ± -> 




1.4a 


B° -» <j)K° s 


o.9±iS 


1.2a 



in the vector meson rest frame fj. 
4.1 B->4>K 

The continuum background suppression is performed by 
sequential cuts; | cos{6 t hr-B)\ < 0.8 and | cos(Q* B )\ < 0.8. 
For B^ — > <^if ± and B° — > channels an additional 
cut is used: |cos(6>#)| > 0.5. The signal yield in each 
studied decay channel is extracted from the unbinned 
maximum-likelihood fit. The fits are made simultane- 
ously in Mt, - AE plane. The signal in M& and AE is 
adequately represented by a single Gaussian. The back- 
ground distribution for Mf, is represented by the ARGUS 
functionf]. The background in AE is parametrized by a 

a In a Pseudoscalar(P) — * Pseudoscalar - Vector(V) decay the vec- 
tor meson is polarized, decaying to two pseudoscalars V —* PP 
according to a cos 2 (0£f); the combinatorial background is flat in 
cos(0h). 

6 B(x) = AVI - a: 2 exp(P(l - x 2 )); where x = E *' b , P is the 

shape variable fitted to background sample and N is the overall 
normalization factor fitted to the distribution in the signal region. 



4.2 B^fm* 

In these channels another discrimination observables are 
used: an angle between the p and the beam axis in the 
B rest frame. In addition the extended Fox- Wolfram] 
moments are used (Super Fox- Wolfram moments ). The 
moments are combined to maximize the separation be- 
tween signal and background distributions. 

The variables are combined into Likelihood ratio C 
for better background rejection. The probabilities of be- 
ing the signal (p sl9 ) and background (p bck ) are calculated 
for each of the variables and the likelihood ratio is formed 
£ = sia i bck - A cut is made at C > 0.9. 

The possible cross talk from B^ — > D°ir ± ;D a — > 
K^ir^ which can imitate B^ — > p a Tr^;p — > n + ir~ is 
removed by checking the D° decay hypothesis. 

The fit for _B ± — > p ^ in the AE projection only 
(with Mb cut) gives the signal yield of 13.73+5;t5 a ^ s i§~ 
nificance of 2.9a and efficiency of 13.4%. The result is 
confirmed in the fit to the Mb distribution, which gives 
a significance of 3.4a. 

For the B° — + / o ± 7r T mode, the cos(6*#) > and 
cos(8h) < regions are fitted separately. The yield is 
determined from a simultaneous fit to the Mb and AE 
distributions 

The fit results are shown in the Tabic ||. The sig- 
nificance of the signal yield predicted by the fit to both 
subsamples combined is 3.5a. 

c The statistical significance of the signal yield is defined as 
\J — 2 \og(Lo / Lmax) , where Lq corresponds to the value at zero 
signal yield hypothesis. 
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Table 4: The fit results for B -* p ± K^ modes. 



data sample 


fitted signal yield 


significance 


cos(9 H ) > 


r n+3.8 
°- U -3.3 


1.9(7 


cos(t9f/) < 


*-°-3.3 


1.5(7 


Table 5: The fit results for B — > u)h modes. 


Mode 


fitted signal yield 


significance 


B ± -> W7T ± 


2.8+j;? 


1.0(7 


B ± -> u/K* 




1.8(7 


B ± -> wh* 


6.4 ±8.9 


1.8(7 



^.3 B -► pi^ ± 

For the B ± — > p°K ± mode the fit is made simultaneously 
in the Mf, and Ai? planes. The fit results are given in 
the Table |[ The significance of the signal is only 2.2<7. 

For the B° — > p ± K T mode the fit in the AE pro- 
jection is performed (with Mf, cut). The fit gives signal 
yield of 8.29+443 with significance of 2. Oct (in the case of 
the fit to Mt„ only the significance is 1.8ct) and efficiency 
of 12.8%. 



Table 6: The summary of the obtained branching fractions and 
upper limits [x!0 — ]. 



Mode 


BF 


UL 90% C.L. 


B+ -> 4>K + 


1 39+0.37+0 14 
1 - oy -0. 33-0.24 


- 


B° -> <j>K° 


- 


< 1.6 


B° -> 


1.5^ ±0.3 


< 3.0 


B+ -> 4>K*+ 


- 


< 3.6 


B+ -» p%+ 


1.12±g;g±0.19 


< 2.88 


B+ p°A'+ 




< 1.35 


5° p+TT- 


2.02+°;^ ± 0.33 


< 3.57 


£° -> p+Jf- 




< 2.36 


5+ -> tuh+ 




< 1.41 


B+ -> 




< 0.94 


5+ -> uoK + 




< 1.05 
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As in the B — ► ph search, a cut on a likelihood ratio 
which includes the uj-amplitudeis used to suppress the 
continuum background. The cut is set at C > 0.95. The 
unbinned likelihood fit to Mb and AE is used to extract 
the signal yields. No signal is observed. 

The urn and u>K yields are obtained from the simul- 
taneous fit. For the u>h mode, the particle identification 
information is not applied in the fast charged tracks se- 
lection. The results are summarized in Table |^. 

5 Conclusions 

Using the data sample of 10.5fb _1 collected by the Belle 
detector, a search was performed for two-body charm- 
less decays including vector meson. A clear signal sig- 
nature of B ± — ► is seen with 7.5ct significance. 
The branching fraction for this mode is measured to be 
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,9+71-+ and B ± -> cfrK** is reported. For B° -> 4>K° S , 
B° -> 4>K*°, B ± -> p ^, B° -> B ± -> LJ7T+ 
and B ± -> wK^ no significant signals are observed and 
the 90% C.L. limits are given. 
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